Method for the optical characterization of materials without using a physical model

ABSTRACT

A method of optical characterization of a layer of material in which the spectrum of index n* (λ) is characterized by a limited number of “nodes” that are points with coordinates (λ i , n i , k i ) or (λ i , n* i ), with n i =n(λ i ), k i =k(λ i ) and n* i =n i +jk i , where j 2   =−1 , and an interpolation law between the “nodes,” which can be, for example, linear, cubic, of “spline” type or polynomial (of any given degree). This interpolation law allows the calculation, from the “nodes,” of the refraction indexes and the extinction coefficients for the wavelengths located between the “nodes.”

TECHNICAL FIELD

The invention relates to a method of optical characterisation of materials.

This method allows to characterise thin or thick layers of these materials, which are created on substrates. The physical dimensions, that this method allows to determine, are:

-   -   the thickness of a layer of material;     -   the refraction index of this material; and     -   the absorption coefficient of this material.

Optical characterisation of materials is used for the chemical analysis of these materials (notably the study of absorption bands, of densification properties and of oxidation properties), in the fields of microelectronics, sensors, biology, medicine, or to analyse the thickness of deposits of these materials.

We refer to document [1] for examples of applications which, as for the other documents cited later on, is mentioned at the end of this description.

The characterisation of the optical properties of a material is also useful when the material is structured later on (to create for example etchings or surface roughness) and the optical diffraction properties of the obtained structure must be calculated (see document [2]).

It can already be noted that the invention is particularly useful when the physical law applied by the complex refraction index of the material to be characterised is, a priori, unknown.

STATE OF THE PRIOR ART

Remember that the optical measurements can have a variety of natures.

It can be reflectometric measurements. In this case, the intensity reflection coefficient of a structure is measured over a spectrum (i.e. an interval) of wavelengths [λ_(m), λ_(M)].

The incident angle of the illuminating light may be not zero. The reflection coefficient can be measured for several incidence angles θ. Let R (θ, λ, p) be the reflectometric spectrum, where p is the polarisation of the incident beam and λ the wavelength of the latter.

Generally, the angle θ is zero and the polarisation p unknown. In the event where θ is not zero, this polarisation p must be known. Generally, the latter is of (S) or (P) type.

It can also be ellipsometric measurements. The measured dimensions are then the real and imaginary parts of the ratio between the polarisation reflection coefficient (P) and the polarisation reflection coefficient (S).

We generally mark ρ=|ρ|exp(jΔ) this complex ratio (with j²=−1). The dimensions commonly handled are |ρ|, that we mark tan(Ψ), and cos(Δ), or combinations of both.

For example, the variables derived from a phase modulation ellipsometer are the following: I _(s)=sin(2Ψ)sin(Δ) and I _(c)=cos(2Ψ).

As for a standard ellipsometer, it supplies the following variables: α=(tan²Ψ−1)/(tan²Ψ+1) and β=cosΔ(1−α²)^(1/2).

For generality purposes, we mark the handled dimensions S₁ and S₂.

The spectrums S_(i), iε [1,2], are measured over a range of wavelengths [λ_(m), λ_(M)]. The incidence angle can be of any value. Several spectrums can be measured at different incidence angles in order to obtain a fuller spectrum. We mark the ellipsometric spectrum s(θ,λ)={S₁(θ,λ),S₂(θ,λ)}.

Additionally, goniometric measurements (reflection coefficient as a function of the incident angle) can be added to the measurement used for characterisation, so as to determine the thickness of the various layers, for one or several wavelengths. These measurements alone are insufficient as we want to determine the complex refraction index over a spectral range from λ_(m) to λ_(M).

So as to simplify the presentation, we mark Ψ a set of reflectometric and/or ellipsometric spectrums (and possibly goniometric spectrums for a few wavelengths).

Without omitting any generality, we describe, in this description, the operating mode of the methods of the prior art and of the invention only for the case of a single thin layer of material, created on a known substrate.

The thickness of this layer is marked ε and the complex refraction index of the material with the wavelength λ is marked n*(λ).

In this regard, remember that the real (respectively imaginary) part of this complex refraction index is marked n(λ) (respectively k(λ)) and called “refraction index” (respectively extinction coefficient”).

Additionally, we mark an error function Er(Ψ⁽¹⁾, Ψ⁽²⁾)—for example the average quadratic gap—between two spectrums Ψ⁽¹⁾ and Ψ⁽²⁾.

For example, we can take, when there are ellipsometric spectrums over several angles θ_(i), iε {1 . . . n}, and a reflectometric spectrum:

$\begin{matrix} {\left. {{{Er}\left( {\Psi^{(1)},\Psi^{(2)}} \right)} = {{\frac{1}{{\,^{\lambda}M} - {\,^{\lambda}m}}{\int_{\lambda\; m}^{\lambda\; M}\left\lbrack {{\frac{1}{n}{\sum\limits_{i = 1}^{n}\;\left\lbrack {{S_{1}^{(1)}\left( {\theta_{i},\lambda} \right)} - {S_{1}^{(2)}\left( {\theta_{i},\lambda} \right)}} \right\rbrack^{2}}} + \left\lbrack {{S_{2}^{(1)}\left( {\theta_{i},\lambda} \right)} - {S_{2}^{(2)}\left( {\theta_{i},\lambda} \right)}} \right\rbrack^{2}} \right\rbrack}} + \left\lbrack {{R^{(1)}(\lambda)} - {R^{(2)}(\lambda)}} \right\rbrack^{2}}} \right\rbrack\ {\mathbb{d}\lambda}} & (1) \end{matrix}$

-   -   with Ψ⁽¹⁾(λ)={S₁ ⁽¹⁾(θ_(i),λ),S₂ ⁽¹⁾(θ_(i),λ),R⁽¹⁾(λ)}     -   Ψ⁽²⁾(λ)={S₁ ⁽²⁾(θ_(i),λ),S₂ ⁽²⁾(θ_(i),λ),R⁽²⁾(λ)}     -   and iε{1 . . . n}

Weighing factors can be introduced into the integral so that the error function can take into account the variations in the accuracy of the measurements of the spectrums.

The optical characterisation of layers of material is generally based around two applications.

The first application is the dimension control of the deposit of thin layers that are used in microelectronics.

Usually we know the deposited material, i.e. we know the complex refraction index of this material at the first wavelengths used for the characterisation.

The laws applied by the complex refraction index are either tabulated or approximated by known physical laws such as, for example, the Cauchy model, the Sellmeier model (see document [3]), the Forouhi laws (see document [4]), and the harmonic oscillator laws (see document [5]). These laws are defined by a finite number of parameters.

For example, a Cauchy type law without absorption, is with two parameters, is defined in the following manner:

${{Re}\left\lbrack {n*(\lambda)} \right\rbrack} = {{n(\lambda)} = {a_{0} + \frac{a_{1}}{\lambda^{2}}}}$ Im[n * (λ)] = k(λ) = 0

When we are certain of the value of the coefficients a_(i)(_(i)ε{0,1}) but we do not know the thickness, a search algorithm is used in order to find the thickness which minimises the error between the measurement Ψ and the theoretical result Ψ taking into consideration the modelled index.

The search algorithm can be, for example, the Simplex method, the Tabou search, the Levendt-Marquart method or the simulated annealing method (see chapter 10 in document [6]).

When the refraction index is approximate, the coefficients a_(i) are integrated into the adjustment process of Ψ and Ψ. The search for the coefficients a_(i) constitutes a method of characterising the refraction index.

However, when the law applied by this refraction index is unknown (it happens that the material is unknown or that it is poorly described by known physical laws), this method remains approximate and the thickness risks being erroneous.

The second application is the characterisation of materials.

The employed method remains unchanged, except that the material is not really known. It is precisely the complex refraction index function that is the closest to reality that is targeted.

The type of law can be chosen via analogy with other materials. However, the law applied by the complex refraction index can be complicated, which is for example the case of a harmonic oscillator law.

$\left\lbrack {{n(E)} + {j\;{k(E)}}} \right\rbrack^{2} = {1 + {\sum\limits_{i = 1}^{n}\;\frac{A_{i}}{E + E_{i} + {j\;\Gamma_{i}}}} - \frac{A_{i}}{E - E_{i} + {j\;\Gamma_{i}}}}$

In the above expression, j²=−1 and the refraction index and the extinction coefficient are expressed as a function of E and not of λ, with E=1240/λ(λ in nm).

In this case, the coefficients of the oscillators are difficult to find if we do not know their value range. The search is difficult to automate, as the search algorithms can give erroneous results and the time wasted can be substantial.

There is an alternative to the search for coefficients: the point to point method (PAP). This PAP method proposes not to choose a physical law and to search the complex refraction index of the material for each wavelength λ_(i), where iε [1 . . . n], with λ₁=λ_(m) and λ_(n)=λ_(M).

For each λ_(i), a search algorithm tries to find the thickness, the index n(λ_(i)) and the extinction coefficient k(λ_(i)) which minimise the error between the measured values Ψ(λ_(i)) and the theoretical result Ψ (λ_(i),n(λ_(i)),k(λ_(i)),ε).

Such a method poses a problem because the different points (λ_(i),ε,n(λ_(i)),k(λ_(i))) are not necessarily physically compatible among themselves: for example, the found thickness can vary depending on the wavelength and the law applied by the complex refraction index, more simply called the index law, can have discontinuities.

This method is usually only valid when the thickness is fully known and that the measurements are of good quality.

PRESENTATION OF THE INVENTION

The purpose of the invention is to resolve the previous inconveniences.

The method which is the object of the invention allows to characterise a material without using a physical model, i.e. without using a physical law applied by the complex refraction index of the material under scrutiny. It is therefore especially useful when such a law is unknown.

This method constitutes an alternative to the aforementioned known methods of characterisation. It can be called “method of nodes” as it uses “nodes”, i.e. points having coordinates (λ_(i), n*_(i)), where n*_(i) is the value taken by the complex refraction index with the wavelength λ_(i) and i takes a limited number of values (whole numbers).

Precisely, the purpose of the invention is a method for optical characterisation of at least one layer of material in an interval A of values taken by a function α of an optical wavelength λ, when λ varies in an interval of wavelengths, this layer being created on a substrate, this method being characterised in that it comprises the following stages:

1) we carry out a group of reflectometry and/or ellipsometry measurements over the interval A, this set of measurements leading to a measured spectrum, marked Ψ, and we choose the methods for calculating associated with the nature of the measurements and with the type of layer to be characterised;

2) we choose m initial values α₁ . . . α_(m) of the function α, belonging to this interval A, m being a whole number at least equal to 1, and we define an interval B as being the set of points α of the interval ranging from the smallest to the biggest number among α₁ . . . α_(m), when m is greater than 1, and as being the interval A when m equals 1;

3) we choose m complex initial values of a complex refraction index n*=n+jk for the m points α_(i), i ranging from 1 to m;

4) when m is not 1, we choose an interpolation law which allows to calculate the refraction index n(α) of the material over the interval B, from the points (α_(i), n_(i)), with n_(i)=n(α_(i)), i ranging from 1 to m, and when m equals 1, n(α) is taken equal to the number n₁(α₁) over the entire interval B;

5) we choose M variable parameters, M being less than or equal to 2 m+1;

6) we choose an error function Er(Ψ, Ψ) which characterises the difference between a measured spectrum Ψ and a theoretical spectrum Ψ;

7) using a minimizing function of Er(Ψ, Ψ) with M parameters, we perform the following series of stages:

a) by applying the interpolation law of (α_(i),n_(i)) is over the interval B, we deduce n(α), α belonging to B;

b) by using n(α) and the thickness ε of the layer, and methods for calculating spectrums, we calculate a theoretical spectrum Ψ (n(α),ε);

c) we compare Ψ and Ψ by using Er(Ψ, Ψ) and, if Er(Ψ, Ψ) is sufficiently small, i.e. less than a predetermined value e, or is minimal, we go to stage e), otherwise we go to stage d);

d) we make the M variable parameters vary so as to tend to the minimum of Er (Ψ, Ψ), and we return to stage a);

e) if Er(Ψ, Ψ) is less than e, we then obtain a set of M variable parameters, for which Er(Ψ, Ψ(n(α,M),ε)) is minimal and the refraction index is then taken equal to the last one obtained, and if Er (Ψ, Ψ) is greater or equal to e we go to stage 8).

8) we increase the number m of initial values of the function α and we return to stage 2).

It is therefore, for example, possible to perform optical characterisation:

-   -   over an interval of wavelengths λ, in this case the interval         [λmin, λmax],     -   or over an interval of inverse wavelengths 1/λ, in this case         over an interval [(1/λ)min, (1/λ)max], where (1/λ)min is equal         to 1/(λmax) and (1/λ)max is equal to 1/(λmin);     -   or over an interval of energies E (with E=hv=hc/λ where h is the         Planck's constant, c the speed of the light in the vacuum and v         the frequency corresponding to (λ), in this case over an         interval [E min, E max], where E min is equal to hc/(λmax) and E         max is equal to hc/(λmin);     -   or, more generally, over an interval [αmin, αmax] of values         taken by a function α of the variable λ.

Furthermore, it is appropriate to note that the invention can be used to characterise a spectrum or a part of a spectrum.

Each interpolation law can be chosen from among the linear interpolation laws, the cubic interpolation laws, the polynomial interpolation laws and the interpolation laws for example of “spline” function type.

According to a preferred embodiment of the method which is the object of the invention, the spectrum is evenly sampled in α (λ), meaning that the initial values of the function α (see the aforementioned stage 2) are evenly distributed over the interval A, the distribution of the nodes thus being homogenous.

As previously seen, α (λ) can be chosen from among λ, 1/λ and hc/λ or any other function of λ, where h is the Plank's constant and c is the speed of the light in the vacuum.

Preferably, in the aforementioned stage 6), we measure the error over an interest interval C which is included in the interval B or equal to this interval B.

The M variable parameters can be the real parts of the refraction indexes at points α_(i), i ranging from 1 to m, or the imaginary parts of these refraction indexes, or these M variable parameters can be constituted by the thickness of the material for which we are searching the refraction index.

Another purpose of the invention is another method for optical characterisation of at least one layer of a material in an interval of wavelengths [λ min, λ max], this layer being created on a substrate, this other method being characterised in that:

-   -   we carry out a set of reflectometry and/or ellipsometry         measurements, this set of measurements leading to a measured         spectrum, marked Ψ;     -   we choose m initial wavelengths λ₁ . . . λ_(m) belonging to this         interval, m being a whole number at least equal to 1, we         associate a refraction index to each wavelength;     -   we choose an interpolation law at least for the refraction index         of the material, for wavelengths lying between the initial         wavelengths λ₁ . . . . λ_(m);     -   we choose M initial parameters, M being at least equal to m,         namaly an initial refraction index n_(i) for each initial         wavelength λ_(i), 1≦i≦m, the initial wavelengths being chosen so         as to determine via interpolation at least the refraction index         for any wavelength within the interval [λ min, λ max], the         couples (λ_(i), n_(i)) being called nodes;     -   we choose reflectometry and ellipsometry methods of calculation;     -   we also choose an error function Er, representative of the         difference between two spectrums Ψ₁ and Ψ₂, the spectrums Ψ₁ and         Ψ₂ being calculated or measured over a number of points greater         than the number m of nodes;     -   using the m initial wavelengths, the M initial parameters and         the interpolation law, we implement the following optimisation         process:     -   we determine a theoretical spectrum, marked Ψ, depending on the         chosen methods of calculation, and on the index deduced via         interpolation of its value at λ_(i), i ranging from 1 to m, over         the spectrum [λmin, λmax];     -   we determine the error Er(Ψ, Ψ), between the measured spectrum         and the theoretical spectrum;     -   we minimise this error by varying the position of the values of         the unknown indexes and/or the thickness of the layer and/or the         values of the refraction indexes with initial wavelengths, and         we obtain a spectrum;     -   we add other wavelengths to the initial wavelengths λi . . . λm,         the added wavelengths constituting new nodes;     -   we repeat the method by choosing a number m′ of initial         wavelengths, m′ being greater than m, and M′ initial parameters,         M′ being greater than M, until the accuracy of each spectrum         thus best represented is equal to a predetermined accuracy.

In this case, according to a first specific embodiment, m is at least equal to 2; according to a second specific embodiment, m is equal to 1 and we choose equal initial refraction indexes.

Also in this case, according to a specific embodiment, the material is non absorbent and the number M is equal to m, the extinction coefficient of the material being set equal to 0; according to another specific embodiment, M is at least equal to 2 m, we furthermore choose an interpolation law for the extinction coefficient of the material, for each initial wavelength λi, 1≦i≦m, furthermore we choose an initial extinction coefficient k_(i), the initial wavelengths furthermore being chosen so as to be able to determine via interpolation the extinction coefficient for any wavelength of the interval [λ min, λ max], and within the optimisation process, we minimise the error by also varying the values of the extinction coefficients at the initial wavelengths and the added wavelengths are furthermore placed so as to best represent the spectrum of the extinction coefficient of the material.

In the case of this other specific embodiment, m can be equal to 1 and we can choose equal initial refraction indexes and equal initial extinction coefficients.

Still in the case of this other method which is an object of the invention, the layer of material can be thin, i.e. with a thickness less than the coherence length of the light used for measuring, we can choose an additional initial parameter, namely an initial layer thickness, and in the optimisation process we can minimise the error by also varying the value of the layer thickness; in an alternative, the layer of material can be thick, i.e. not thin, and M can be at most equal to 2 m; in another alternative, the thickness of the layer of material can be known with sufficient accuracy and M is at most equal to 2 m.

The distribution of the nodes can be homogenous.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention will be better understood upon reading the description of the below embodiments, given by way of non-restrictive examples, making reference to the annexed drawings among which:

FIG. 1 is a diagrammatic view of devices allowing to characterise a layer according to the invention;

FIG. 2 shows the variations in the refraction index as a function of the wavelength, for a material obeying to a Cauchy law (curve I) and for a material characterised according to the invention (curve II);

FIGS. 3A and 3B show the variations in the refraction index (respectively in the extinction coefficient) as a function of the wavelength, for a material obeying to a law with two harmonic oscillators (curve I) and for a material characterised according to the invention; and

FIG. 4 diagrammatically illustrates the parameters used in a generalisation of examples of the invention.

DETAILED DESCRIPTION OF PARTICULAR EMBODIMENTS

The invention proposes an alternative to the aforementioned standard methods. It allows to combine the coherence of a layer model (corresponding to a continuous index law and to constant physical thicknesses), with the generality concerning the index law to be found (as in the PAP method). Furthermore, the resolution is only limited by the resolution of the measured spectrum.

In the method according to the invention, the spectrum of index n*(λ) is characterised by:

-   -   a limited number of “nodes” which are points having co-ordinates         (λ_(i), n_(i), k_(i)) or (λ_(i), n*_(i)), with n_(i)=n (λ_(i)),         k_(i)=k (λ_(i)) and n*_(i)=n_(i)+jk_(i), where j²=−1, and     -   an interpolation law between the nodes, which can be, for         example, linear, cubic, of “spline” type or polynomial (of any         given degree).

This interpolation law allows to calculate, from the nodes, the refraction indexes and the extinction coefficients for the wavelengths located between the nodes.

For example, when the refraction index is characterised by a set of value for wavelengths λ₁ . . . λ_(m), we can use a linear interpolation between two wavelengths λ_(i) and λ_(i+1) to calculate the index n at the wavelength λ (see document [6] chapter 3):

${n(\lambda)} = \frac{{\left( {\lambda_{i + 1} - \lambda} \right){n\left( \lambda_{i} \right)}} + {\left( {\lambda - \lambda_{i}} \right){n\left( \lambda_{i + 1} \right)}}}{\lambda_{i + 1} - \lambda_{i}}$

-   -   with λ_(i)<λ<λ_(i+1)

We can do the same for the extinction coefficient.

When the number of nodes makes it possible, more complex interpolation formulae, bringing into play nearby nodes, can be used (see document [6] chapter 3).

A layer model is therefore characterised by a thickness ε and a family of nodes.

Hereafter we present an example of the method according to the invention. In this example, the measurements Ψ are constituted of a reflectometry measurement R(λ) and an ellipsometry measurement S_(1,2)(θ,λ) where θ is the incidence angle of the light beam that is transmitted onto the layer to be studied during measuring via ellipsometry.

This layer is a thin layer so that the thickness of this layer is also a variable of the problem. Moreover, we suppose that only one layer is unknown, this layer being created on a known substrate.

We will first of all briefly explain this example which uses an algorithm (algorithm of the “nodes method” according to the invention).

From assumed information about the thickness ε of the scrutinised layer and about the refraction index n(λ) and about the extinction coefficient k(λ) of the material of this layer, we create initial nodes (only a few) and an initial thickness ε.

We thus have m nodes and, via interpolation, we can find n(λ) and k(λ) when λ is different from the values of the wavelengths associated with the nodes.

From the initial thickness ε and these initial values n(λ) and k(λ), we determine the theoretical spectrum Ψ by using ellipsometric and reflectometric calculations.

Furthermore, by means of ellipsometric and reflectometric devices and a spectrometer, we obtain S_(1,2)(θ,λ) and R(λ) and we deduce the measurements marked Ψ (for measuring conditions θ and λ).

We then compare Ψ and Ψ by using an error function Er and we optimise the refraction index value and the extinction coefficient value at the different nodes as well as the thickness value, by trying to minimise Er(Ψ, Ψ).

Once these values have been optimised and if the accuracy on the spectrum n(λ), the spectrum k(λ) and the thickness ε is insufficient, we introduce new nodes, we make the thickness ε vary and we recommence the determination of Ψ, the comparison of Ψ and Ψ and the optimisation that are aforementioned, etc.

We break the thus defined loop when the accuracy on each of the spectrums n(λ) and k(λ) and on the thickness ε is considered to be sufficient (satisfactory adjustment of Ψ and Ψ).

The spectrums n(λ), k(λ) and the thickness ε are thus characterised.

In a non-restrictive illustration we find a layer thickness ε equal to 212.3 nm.

FIG. 1 diagrammatically shows the scrutinised layer 2, created on a substrate 4. We can see the ellipsometric device 5,6, the reflectometric device 8 and the spectrometer 10. Furthermore, we can see electronic processing means 12, comprising a computer and allowing to characterise n(λ), k(λ) and ε according to the information supplied by the spectrometer 10 and according to the method of the invention.

These means 12 are equipped with display means 14 which allow, in particular, to display the curve of the variation of n as a function of λ and the curve of the variation of k as a function of λ.

We will now consider the illustration in greater detail.

Phase 1

The method of this example firstly comprises an initialisation stage.

The algorithm starts with a limited number of nodes, more precisely at least one node. We can therefore start with a single node, by imposing a refraction index and an extinction coefficient which remain constant when the wavelength varies.

We chose node positions so as to, from this node family, deduce the entire spectrum via interpolation. The layer model therefore has 3 parameters or more, as the thickness is also a variable to be determined. The index table over the entire spectrum is therefore deduced from the nodes via interpolation.

This is the case when the thickness of the layers is about the dimension of the wavelength (thin layers). But when the thickness of the layer is greater than the coherence length of the light source, the thickness barely enters into the calculation of the layer response, and is therefore no longer a variable of the problem.

For example, in the case of optical disks (CDROM), on which a very thick layer is deposited (a thickness of about 1 millimetre), the reflection coefficient of such a layer is no longer a function of the thickness of the layer but purely of the refraction index of the latter, the coherence length of the incident light beam being less than the thickness of the layer. In this precise case, the coherence length of the incident beam is determined by the roughness of the layers.

We choose, for example, to place the first two nodes at the ends λ min and λ max of the spectrum. The complex index values at these ends are chosen according to the type of material under scrutiny. For example, on an ellipsometric spectrum between 300 nm and 800 nm of a thin layer of photoresist, we take n(300 nm)=n(800 nm)=1.5 and k(300 nm)=k(800 nm)=0.

When the spectrum is only characterised by two nodes, the index between the ends is determined via linear interpolation. In the case under consideration, we therefore have n(λ)=1.5 and k(λ)=0 for belonging to [300 nm, 800 nm].

From three nodes up, we rather opt for cubic interpolation in order to obtain softer forms of index law compared to those obtained via linear interpolation.

The initial thickness is, itself, chosen as close as possible to the real thickness.

Phase 2

We then proceed with an optimal determination of the values of the refraction index and of the extinction coefficient on the nodes and of the thickness value.

To accomplish this, the spectrums Ψ(λ) are calculated by means of the used layer model, resulting from the choice of nodes, the interpolation law and the thickness of the layer.

The physical model used to calculate Ψ of course depends on the applied method for measuring, i.e. in particular on the incidence angle of the light, on the applied spectrum and on the model of thin layers or of thick layers if appropriate (see for example the model of stacked layers in document 3)).

The spectrums Ψ being constituted of a set of measurements of diverse nature (for example ellipsometric and reflectometric measurements), we use a reflectometric (respectively ellipsometric) calculation method for reflectometric (respectively ellipsometric) measurements.

The reflectometric and ellipsometric measurements are combined by means of an error function Er(Ψ, Ψ) which is for example of the type defined by the equation (1).

Thanks to a search function, we minimise the difference between Ψ(λ) and Ψ(λ), by varying the value of the refraction index and the value of the extinction coefficient at the position of each of the nodes as well as the layer thickness (if this thickness is an important factor in the calculation of Ψ). When the difference is minimal, that being when Er(Ψ, Ψ) is minimal, this signifies that the reflectometric and ellipsometric measurements coincide to their best (for a given number of nodes).

At this stage we obtain, for a known number of nodes and a known spectral position for each of these nodes, the layer model (refraction index, extinction coefficient and thickness) which best corresponds to the real layer.

The validity of the found model is better ensured as the number of measurements increases. To obtain a large number of measurements we can for example use several incidence angles of the light θ_(i), 1≦i≦l, for the ellipsometric spectrums, carry out a reflectometric measurement and complementary goniometric measurements.

Phase 3

Then, we increase the number of nodes.

We add a finite number of nodes. In a first embodiment, the added nodes are positioned so as to best represent the spectrums n(λ) and k(λ). By way of example, we place these additional nodes at spots where the difference between Ψ and Ψ is maximum or at spots where the nodes are the furthest apart. We then return to phase 2 as long as the accuracy on each of the spectrums n(λ) and k(λ) and on the thickness ε is insufficient, i.e. not equal to a pretedermined accuracy.

In a second embodiment, we insert new nodes between two nodes of the set of previously chosen nodes, these new nodes being evenly distributed over the spectrum.

All the same it is appropriate to heed the following.

When we increase the number of nodes, the position of the old nodes does not need to be conserved. For example, assume that we evenly sample a spectrum from 400 nm to 800 nm with 3 nodes. These nodes are therefore respectively located at 400 nm, 600 nm and 800 nm. When we increase to 6 nodes, the three additional nodes cannot be placed so that the spectrum is evenly sampled if we conserve the position of the old nodes. The position of the 6 nodes can be defined, is we wish to have even sampling, by the values 400, 480, 560, 640, 720 and 800 nm. The old middle node at 600 nm therefore disappears. To calculate the index value at these positions from the old nodes we use interpolation.

In the following we provide two common examples of application of the invention. These examples implement two types of materials which apply different laws.

From ellipsometric and reflectometric measurements we propose to find the physical laws applied by these materials.

We proceed in the following manner.

A fictitious material is created, this material applying a known theoretical law (a Cauchy law or a harmonic oscillator law), with parameters that we arbitrarily set. The variations of the complex refraction index as a function of the wavelength are thus perfectly known.

Furthermore, we impose a thickness of 200.00 nm for the material, on a silicon substrate, this latter also being well known.

Fictitious measurements (ellipsometric measurements, reflectometric measurements) are calculated, then noise induced so as to introduce an apparatus defect.

It all takes place as if we had used real measurements taken on the material. But contrary to reality we fully know the complex refraction index as we set it just as we set the thickness of the layer of material.

Here we “blindly” test the method, that meaning we start with an erroneous thickness (220 nm) and false complex refraction indexes, as they are supposed to be unknown.

We apply the method according to the invention then we compare the complex refraction index found with the is theoretical complex refraction index. We precisely find the same laws, as well as the same layer thickness.

As a first example, take a material whose complex refraction index applies to a Cauchy law such as:

${n(\lambda)} = {1.5 + {0.1\frac{300^{2}}{\lambda^{2}}} + {0.1\frac{300^{4}}{\lambda^{4}}}}$ k(λ) = 0

This index law is typical of photoresists (in the spectral range from 300 nm to 800 nm).

In order to find, using the nodes method (i.e. the method according to the invention), the aforementioned index law, we carry out two measurements, namely an ellipsometric measurement at a 70° angle and a reflectometric measurement.

The processing conditions are as follows:

-   -   the processed spectrum lies between 300 nm and 800 nm;     -   at the start, the nodes are in the positions (300 nm, 1.6) and         (800 nm, 1.6), this meaning that the index is considered as         linearly varying between 300 nm and 800 nm, and that its value         is constant (equal to 1.6);     -   the number of nodes is repeatedly increased according to the         sequence 2→4→6;     -   during the node increasing procedure, the wavelength position of         each node is calculated so that the sampling by 1/λ is regular         (λ: wavelength), the density of points thus being increased         around the low wavelengths;     -   the interpolation law is a cubic law, when the number of nodes         is greater than 2, otherwise it is linear; and     -   the applied minimisation algorithm is a Simplex type algorithm.

FIG. 2 allows to compare the refraction index corresponding to the fictitious material which perfectly applies the Cauchy law (curve I) with the refraction index which we find via the nodes method (curve II), using 6 nodes (represented by circles in FIG. 2). We have used a Ψ composed of an ellipsometric measurement {S₁(λ),S₂(λ)} at 70° and a reflectometric measurement R(λ).

We finally find a thickness of 199.8384 nm.

As for a second example, take a material whose complex refraction index applies a law with two harmonic oscillators, such as:

$\left\lbrack {{n(E)} + {j\;{k(E)}}} \right\rbrack^{2} = {1 + {\sum\limits_{i = 1}^{2}\;\frac{A_{i}}{E + E_{i} + {j\; G_{i}}}} - \frac{A_{i}}{E - E_{i} + {j\; G_{i}}}}$

with

j²=−1 and E=1240/λ (λ in nm)

A₁=0.25×1.5² A₂=0.25×1.5²

E₁=1240/400 E₂=1240/300

G₁=0.3 G₂=0.3

In this second example, the nodes method is applied to a set of ellipsometric measurements carried out between 250 nm and 800 nm, at 75°, 70°, 60° and 45°, along with a reflectometric measurement. The real thickness of the material being 200 nm, we find a thickness of 200.25 nm with the nodes method. The adjustment to the index law considered in this second example is very good, as FIGS. 3A and 3B demonstrate.

These FIGS. 3A and 3B respectively illustrate the reconstructions of the curves n(λ) and k(λ) of the material by means of the nodes method. The reconstruction is carried out using four ellipsometric spectrums and a reflectometric spectrum. The real absorption peaks are particularly well represented by the curve obtained via cubic interpolation between the nodes (represented by the circles in FIGS. 3A and 3B).

In FIG. 3A, the curve I (respectively II) corresponds to a refraction index n which perfectly applies the chosen law (respectively to a refraction index n found via the nodes method).

In FIG. 3B, the curve I (respectively II) corresponds to an extinction coefficient k which perfectly applies the chosen law (respectively to an extinction coefficient k found via the nodes method).

Examples of the invention have just been described above. More generally speaking, we note that in the invention we consider a set X of values, with X={n₁,n₂, . . . ,n_(i), . . . n_(m),k₁,k₂, . . . k_(i), . . . k_(m),ε}, where

n_(i) is the value of the refraction index (real) corresponding to λ_(i), iε{1 . . . m}, m being the number of nodes

k_(i) is the value of the absorption coefficient at the node corresponding to λ_(i), iε{1 . . . m},

ε is the thickness of the scrutinised layer.

In this case, the operation for minimizing the error Er(Ψ, Ψ) means finding the set or “vector” X such that Er is minimal.

When we do not assert any specific constraints, the minimizing is a minimizing with 2×m+1 parameters. We can of course introduce constraints in order to reduce the number of variables.

In particular, if we know that the material is non-absorbent, we impose k_(i)=0 for all i of {1 . . . m} and X becomes: X={n₁,n₂ . . . ,n_(i), . . . n_(m),ε}.

If through a complementary measurement (for example a goniometric measurement or a non-optical direct measurement) we know with sufficient accuracy the thickness of the scrutinised layer, the thickness ε is no longer a variable and we obtain: X={n_(i), n₂ . . . , n_(i), . . . n_(m), k₁, k₂, . . . k_(i), . . . , k_(m)}.

Of course, the two previous options can be combined.

In the following we explain an embodiment of the invention in a more general form than the previous examples.

Given that α(λ) is a function of the wavelength λ of the light used in the measurements, we can for example choose:

α(λ)=(see FIGS. 3A and 3B where the spectrum is evenly sampled in λ)

α(λ)=1/λ (see FIG. 2 where the spectrum is evenly samples in 1/λ)

α(λ)=hc/λ where h is the Planck's constant and c the speed of light in the vacuum, α(λ) then being homogenous with an energy.

Given that A is the interval of the spectrum of measurement, B is the spectral interval described by the “nodes” and C is the interval of interest.

The interval C is included in the interval B or equal to this interval B. Likewise, the interval B is included in the interval A or equal to this interval A.

We specify that each of the intervals A, B and C is is of [α_(m), α_(M)] type where α_(m) is less than α_(M) and there are two wavelengths λ_(k) and λ_(i) such that α_(m)=α(λ_(k)) and α_(M)=α(λ₁). By way of non-restrictive illustration, FIG. 4 shows an example of the intervals A, B and C and of the variation curves of Ψ and Ψ as a function of α(λ), Ψ in fact being a function of n*(α(λ)).

The circles N represent the nodes.

n* is the complex index which is expressed here as a function of α(λ) and whose real and imaginary parts are respectively marked n(α(λ)) and k(α(λ)).

We also see, by way of non-restrictive example, an example of the variations curve of n (respectively k) as a function of α(λ), passing through the points having co-ordinates (n_(i),α_(i)) (respectively (k_(i), α_(i)), where α_(i)=α(λ_(i)), n_(i)=n(α_(i)), k_(i)=k(α_(i)), 1≦i≦m (m being a natural number other than zero).

In FIG. 4 we note the relation between the nodes N, the points (n_(i), α_(i)) and the points (k_(i), α_(i)).

In the embodiment under consideration, we use an algorithm comprising the following stages:

1. we carry out the measurements Ψ over the interval A and we choose the calculation methods associated with the measurements (ellipsometric or reflectometric calculations);

2. we choose m numbers α_(i) (constituting m initial values of the function α), i belonging to {1, . . . m}, with m≧1, and {α_(i)}⊂A (the α_(i) corresponding to the “nodes”); when m>1, B is defined as the set of points α such that min(α_(i))≦α≦max(α_(i)); when m=1, we have B=A;

3. we choose m initial values of complex index n* at the m points α_(i), i belonging to {1, . . . m};

4. if m≠1 we choose an interpolation law which allows to calculate the refraction index n(a) over the interval B from the points (α_(i), n_(i)), i belonging to {1, . . . , m}; if m=1, then n(α)=n₁(α₁) over the entire interval B;

5. we choose M variable parameters with M≦2 m+1; these parameters can be for example:

-   -   the real parts of the refraction indexes at points α_(i), i         belonging to {1, . . . , m}, or     -   the imaginary parts of the refraction indexes at these same         points, or     -   the thickness of the material for which we are searching the         refraction index;

6. we choose an error function Er(Ψ, Ψ) which characterises the difference between a measured spectrum and a theoretical spectrum; generally, the error is measured over the interval C;

7. using a minimizing function of Er(Ψ, Ψ) with M parameters, we perform the following series of stages:

a) by applying the interpolation function of (α_(i),n_(i)) over B, we deduce n(α) with α belonging to B;

b) by using n(α) and the thickness ε and methods for calculating spectrums, we calculate the theoretical spectrum Ψ (n(α), ε);

c) we compare Ψ and Ψ by using Er=Er(Ψ, Ψ); if Er is sufficiently small (i.e. if Er is less than a predetermined value e), or if Er is minimal, we go to stage e), otherwise we go to stage d);

d) we make the M variable parameters vary so as to tend to the minimum of Er and we return to stage a);

e) if Er is less than e we then obtain a set of M parameters such that Er(Ψ, Ψ(n(α,M),ε)) is minimal and the refraction index calculation is completed: this index is taken equal to the last one obtained; and if Er is greater than or equal to e we go to stage 8).

8. We increase m and we return to stage 2).

It is appropriate to note that the invention can not only use wavelength (λ) samplings but also frequency (c/λ or simply 1/λ) samplings, energy (hc/λ) samplings and, more generally, samplings with parameters which are functions of the wavelength.

Furthermore, it is appropriate to note that an aforementioned essential stage of the algorithm (stage 8) is not restricted to adding a set of nodes to the already existing nodes: it encompasses the more general case where the number of nodes increases. This means, in a specific embodiment of the invention, that after minimizing with 3 nodes, if we want to increase to 6 nodes in all, the position of the 3 old nodes is “erased” so as, for example, to have a constant density of nodes over the spectrum. In practice, this is the best option. The information about the position of the old nodes is not lost as the index value at the old nodes is used to calculate the values of the 6 new nodes (in fact 3 new nodes plus 3 old nodes).

Thus, according to a specific embodiment of the method according to the invention, we can increase the number of initial values of the function α by adding one or several values to the extant initial values; however, according to a preferred embodiment, we can increase the number of initial values of the function α by replacing the extant initial values with new initial values whose number is greater than the number of extant initial values.

The invention in not restricted to the characterisation of thin layers. It also applies to the characterisation of thick layers.

Furthermore, the invention is not restricted to the characterisation of a single layer, created on a substrate. It also applies to the characterisation of two, or more than two, layers created on a substrate.

The following documents were cited in the description:

[1] R. M. A. Azzam and N. M. Bashara, Ellipsometry and Polarised Light, North-Holland Physics Publishing, 1997, chapter 6.

[2] B. K. Minhas, S. A. Coulombe, S. Sohail, H. Naqvi and J. R. McNeill, Ellipsometry scatterometry for the metrology of sub-0.1-μm-linewidth structures, Applied Optics, 37(22): 5112-5115, 1998.

[3] M. Born and E. Wolf, Principle of Optics, Cambridge University Press Edition.

[4] A. R. Forouhi and I. Bloomer, Optical dispersion relations for amorphous semiconductors and amorphous dielectrics, Physical Review B, 34(10): 7018-7026, November 1986.

[5] F. L. Terry, Jr., A modified harmonic oscillator approximation scheme for the dielectric constants of Al_(x)Ga_(1-x)As, Journal of Applied Physics, 70(1), 1991, pages 409-417.

[6] W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in C, Cambridge University Press, 1992, chapters 3 and 10. 

1. A method for optical characterization of at least one layer of material in an interval A of values taken by a function α of an optical wavelength λ, when λ varies in an interval of wavelengths, the at least one layer being created on a substrate, the method comprising: 1) carrying out a set of reflectometry and/or ellipsometry measurements over the interval A with ellipsometric and/or reflectometric devices and a spectrometer, the set of measurements leading to a measured spectrum, marked ψ, and choosing methods for calculating associated with a nature of the measurements and with a type of layer to be characterized; 2) choosing m initial values α₁ . . . α_(m) of the function α, belonging to the interval A, m being a whole number at least equal to 1, and defining an interval B as being the set of points α of the interval ranging from the smallest to the biggest number among α₁ . . . , α_(m), when m is greater than 1, and as being the interval A when m equals 1; 3) choosing m complex initial values of a complex refraction index n*=n+jk for the m points α_(i), i ranging from 1 to m; 4) when m is not 1, choosing an interpolation law that allows calculating the refraction index n(α) of the material over the interval B, from the points (α_(i), n_(i)), with n_(i)=n(α_(i)), i ranging from 1 to m, and when m equals 1, n(α) is taken equal to the number n1(α1) over the entire interval B; 5) choosing M variable parameters, M being less than or equal to 2 m+1; 6) choosing an error function Er (ψ, ψ) that characterizes the difference between a measured spectrum ψ and a theoretical spectrum ψ; 7) using a minimizing function of Er (ψ, ψ) with M parameter, performing: a) by applying the interpolation law of (α_(i), n_(i)) over the interval B, deducing n(α), α belonging to B; b) by using n(α) and the thickness ε of the layer, and methods for calculating spectrums, calculating a theoretical spectrum ψ(n(α),ε); c) comparing ψ and ψ by using Er(ψ, ψ) and, if Er(ψ, ψ) is less than a predetermined value e, or is minimal, going to sub-step e), otherwise going to sub-step d); d) making the M variable parameters vary so as to tend to the minimum of Er(ψ, ψ), and returning to sub-step a); e) if Er(ψ, ψ) is less than e, then obtaining a set of M variable parameters, for which Er(ψ, ψ(n(α,M),ε)) is minimal and the refraction index is then taken equal to the last one obtained, and if Er(ψ, ψ) is greater or equal to e going to step 8); 8) increasing the number m of initial values of the function α and returning to step 2).
 2. A method according to claim 1, further comprises increasing the number of initial values of the function a by adding one or plural values to extant initial values.
 3. A method according to claim 2, further comprising increasing the number of initial values of the function a by replacing the extant initial values with new initial values whose number is greater than the number of extant initial values.
 4. A method according to claim 2, wherein each interpolation law is chosen from among linear interpolation laws, cubic interpolation laws, polynomial interpolation laws, and interpolation laws of spline function type.
 5. A method according to claim 1, wherein the initial values of the function a are evenly distributed over the interval A, the distribution of the nodes thus being homogenous.
 6. A method according to claim 1, wherein α(λ) is chosen among λ, 1/λ and hc/λ, where h is the Planck's constant and c the speed of light in vacuum.
 7. A method according to claim 1, further comprising measuring the error, at 6), over an interest interval C which is included in the interval B or equal to the interval B.
 8. A method according to claim 1, wherein the M variable parameters are real parts of the refraction indexes at points α_(i), i ranging from 1 to m.
 9. A method according to claim 1, wherein the M variable parameters are imaginary parts of the refraction indexes at points α_(i), i ranging from 1 to m.
 10. A method according to claim 1, wherein the M variable parameters are constituted by the thickness of the material for which the refraction index is searched.
 11. A method for optical characterization of at least one layer of a material in an interval of wavelengths [λ min, λ max], the at least one layer being created on a substrate, the method comprising: carrying out a set of reflectometry and/or ellipsometry measurements with ellipsometric and/or reflectometric devices and a spectrometer, the set of measurements leading to a measured spectrum, marked ψ; choosing m initial wavelengths λ₁ . . . λ_(m) belonging to the interval, m being a whole number at least equal to 1, and associating a refraction index to each wavelength; choosing an interpolation law at least for the refraction index of the material, for wavelengths lying between the initial wavelengths λ₁ . . . λ_(m); choosing M initial parameters, M being at least equal to m, an initial refraction index n_(i), for each initial wavelength λ_(i), 1≦i≦m, the initial wavelengths being chosen so as to determine via interpolation at least the refraction index for any wavelength within the interval [λmin, λmax], couples (λ_(i), n_(i)) being nodes; choosing reflectometry and ellipsometry methods of calculation; choosing an error function Er, representative of the difference between two spectrums ψ₁ and ψ₂, the spectrums ψ₁ and ψ₂ being calculated or measured over a number of points greater than the number m of nodes; using the m initial wavelengths, the M initial parameters, and the interpolation law, implementing an optimization process of: determining a theoretical spectrum, marked ψ, depending on the chosen methods of calculation, and on the index deduced via interpolation of its value at λ_(i), i ranging from 1 to m, over the spectrum [λmin, λmax]; determining the error Er (ψ, ψ), between the measured spectrum and the theoretical spectrum; minimizing the error by varying the position of the values of the unknown indexes and/or the thickness of the layer and/or the values of the refraction indexes with initial wavelengths, and obtaining a spectrum; adding other wavelengths to the initial wavelengths λi . . . λm, the added wavelengths constituting new nodes; repeating the method by choosing a number m′ of initial wavelengths, m′ being greater than m, and M′ initial parameters, M′ being greater than M, until the accuracy of each spectrum thus best represented is equal to a predetermined accuracy.
 12. A method according to claim 11, wherein m is at least equal to
 2. 13. A method according to claim 11, wherein m is at least equal to 1 and equal initial refraction indexes are chosen.
 14. A method according to claim 11, wherein the material is non absorbent and the number M is equal to m, the extinction coefficient of the material being set equal to
 0. 15. A method according to claim 11, wherein: M is at least equal to 2 m; the method further comprising: choosing an interpolation law for the extinction coefficient of the material; each initial wavelength λi, 1≦i≦m choosing an initial extinction coefficient k_(i), the initial wavelengths furthermore being chosen so as to be able to determine via interpolation the extinction coefficient for any wavelength of interval [λ min, λ max]; within the optimization process, minimizing the error by also varying values of the extinction coefficients at the initial wavelengths, and further placing the added wavelengths so as to best represent the spectrum of the extinction coefficient of the material.
 16. A method according to claim 15, wherein m is equal to 1 and equal initial refraction indexes and equal initial extinction coefficients are chosen.
 17. A method according to claim 11, wherein the layer of material has a thickness less than coherence length of light used for measuring, and further comprising choosing an additional initial parameter of an initial layer thickness, and in the optimization process the error is minimized by also varying the value of the layer thickness.
 18. A method according to claim 11, wherein the layer of material is thick, and M is at most equal to 2 m.
 19. A method according to claim 11, wherein the thickness of the layer of material is known with a predetermined accuracy and M is at most equal to 2 m.
 20. A method according to claim 11, wherein each interpolation law is chosen from among linear interpolation laws, cubic interpolation laws, polynomial interpolation laws, and interpolation laws for example of spline function type.
 21. A method according to claim 11, wherein a distribution of the nodes is homogenous. 